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1 Abstract 

In this letter we describe a method of extending an existing phase field model 
of polycrystalline solidification from two to three dimensions (3D). 

2 Introduction and Preliminaries 

There have been a number of approaches to the modeling of grain boundaries, 
all of which have limitations and advantages. Of particular interest are phase 
field models, which have gained popularity as their ability to compute realistic 
microstructures has been demonstrated. For an overview of this approach, the 
reader is recommended some of the review articles on this topic. pQ 

A basic model of grain boundaries in 2D (see can be derived from the 

total free energy 

f a 2 e 2 

T = I dV[f(</>, T) + — \V<j>\ 2 + sg(cf>)\Ve\ + — /i(0)|V^| 2 ] , 

where /(</>, T) + ^-|V</>| 2 are the terms found in classical phase field models of 
solidification, namely the bulk free energy density, which depends on the phase 
field <f> and the temperature T, with minima in the liquid and solid phases 
4> — 0, 1 plus a gradient penalty for interfaces. For this discussion we have 
omitted terms accounting for interface energy anisotropy, although such effects 
are both important and can be accounted for with well known extensions to this 
theory. |2] 

The final terms in the free energy arc functions of the gradient in the orienta- 
tion, 9; introduced to allow for grain boundary energy misoricntation penalties. 
These terms inclusion also provides a realistic description of related phenomena 
such as polycrystalline growth and nucleation.0],|S] The couplings g((j>) and h((j)) 
are chosen so there are no energy penalties in the liquid (i.e. g(0) = h(0) = 0.) 
The dynamics of the system are found by imposing the thermodynamic require- 
ment that <j) and 9 evolve so as to minimize the free energy F. 



1 



3 Extension to 3D: Formulation and Solution 



We note that a single angle cannot represent an orientation in 3D, and thus 
this concept must be replaced with a more robust mathematical description 
of orientation. Specifically, 9 must be replaced with an object that captures 
the three rotational degrees of freedom available in 3D. Additionally, we must 
also define the norm of this object. With these two mathematical concepts the 
transition to three dimensions is fully posed. 

3.1 Formulation 

There are many ways to represent orientations in 3D, most quite familiar to crys- 
tallographers: Euler angles, rotation vectors, Rodrigues vectors, quaternions, 
etc. All of these representations are mathematically equivalent representations 
of the group SO(3) (special orthogonal group 3D), but retain advantages and 
disadvantages, depending on the application. If we call a member of this group 
P, then P is a 3 x 3 orthogonal matrix (P T P = /, where / is the identity matrix) , 
and P has a positive determinant det P = 1. Thus, we say 50(3) is naturally 
embedded in R 9 , as it can be represented as a nine-dimensional object (with 6 
constraints originating from the orthogonality condition). 

To proceed, we must find the 3D analog to the fundamentally 2D quantity 
|V0|. A gradient is simply a difference over an infinitesimal distance, thus we 
need to compute the norm of the difference, in some sense, between two 3D 
orientations. We consider two possible choices for this measure, employing a 
function of two SO(3) matrices p(P,Q): 

Type I : p(P, Q) = \P-Q\ = {PQ^ 1 - I\ 

Type II : p(P, Q) = \/2 cos' 1 tlP Q~^ ~ 1 

Note that, PQ' 1 is the misorientation between two crystals. The meaning of 
the Type I measure is trivial, as it measures the distance between two matrices 
in R 9 . The Type II measure, on the other hand, measures the length of the 
geodesic in 5*0(3) connecting two matrices. These two measures coincide when 
P and Q are infinitesimally close, but they will yield different values when there 
is a discontinuity between P, Q as is often the case for discrete computations on 
a lattice 

With these definitions we can now write down our model, by simply substi- 
tuting |V0| -> \VP\ in Eqn. □ 

T = [ dV[f(<P, T) + ^|V0| 2 + S g(4>)\VP\ + f M0)|Wf] , 



A more explicit form can be obtained using \VP\ = y/\ VP\ 2 — yY^j=i \^Pi,j\ 2 i 
where p it j = [P]i,j. 
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3.2 Solution 



Having posed the above free energy, we must now perform a minimization, 
to derive equations of motion for both the phase field <f> and the orientation. 
We must proceed with care to ensure that the equations of motions keep the 
variables in SO(3) . There are several ways to proceed: (1) Derive equations for 
the constrained free energy, which has only 3 degrees of freedom or (2) derive 
equations on R 9 , and project the results back into SO(3). For both methods 
we need to derive the variational derivative of the free energy with respect to 
orientation, which is simply 

g = -v-(^)^,+^)vp). 

In deriving the equations of motion for the constrained free energy, an ele- 
ment of SO(3) is written in the form P = P(u, v, w) where the triplet (u, v, w) 
is some local coordinate, for example, the Rodrigues vector. The equations of 
motion for these variables should be 

du _ IJ_T_ dP\ 

Tu dt~\ SP> du) { ) 

with identical equations for u — > v, w. Note that the quantity (•, •) is the usual 
inner product in R 9 (a fully contracted matrix product), and t u is an inverse 
mobility. 

Alternatively, to use a projective formulation, we develop 9 equations of 
motion in R 9 keeping the solution within SO(3) by taking a projection of driving 
force onto the tangential plane of SO(3). It is given in the form TpdP/dt = 
np (—SJ-/SP) where np is the projection operator. This approach allows for 
substantially improved numerical efficiency However, we reserve discussion of 
this technique for a later publication. 

Following the preceding arguments we can derive the evolution equation 
for and P, implement the equations in computer code and solve. Our first 
calculations were for a thin film where the grains are nearly 2D objects, but 
their orientation is 3D, and the dynamics will be governed by the evolution 
of all of these angles. There are numerous experimental systems analogous 
to this calculation (see Fig^) Our results showing growth, impingement and 
coarsening are given in Fig. [2J 

Herein, we have extended our previous work in 2D to 3D. For a complete 
formulation to be obtained, we must also include the important consequences 
of anisotropy. In other words we have examined the consequences of misori- 
entation, but not the consequences of inclination on the statics and dynamics 
of grain boundaries. Additionally, we have yet to account for the underlying 
crystal symmetries These effects can be included, and will be discussed in future 
work. 
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Figure 1: Coarsening process of grain structure of succinonitrile. The grain 
structure is almost 2D, while the orientation of each grain is necessarily 3D. 
(courtesy of Drs. Lee and Losert, U. Maryland) 




Figure 2: Simulation of solidification and coarsening process (the color indicates 
one of the three Euler angles (all of which were solved for). 
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